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COMPUTATIONAL  DETERMINATION  OF  HEATS  OF  FORMATION 
OF  ENERGETIC  COMPOUNDS 

PETER  POLITZER,  JANE  S.  MURRAY  AND  M.  EDWARD  GRICE 
University  of  New  Orleans,  Department  of  Chemistry,  New  Orleans,  Louisiana 

70148  USA 

ABSTRACT 

A  recently-developed  densitv  functional  procedure  for  computing  gas  phase 
heats  of  formation  is  briefly  described  and  results  for  several  categories  of  energetic 
compounds  are  summarized  and  discussed.  Liquid  and  solid  phase  values  can  be 
obtained  by  combining  the  gas  phase  data  with  heats  of  vaporization  and  ^ 
sublimation  estimated  by  means  of  other  relationships.  Some  observed  functional 
group  effects  upon  heats  of  formation  are  noted. 

INTRODUCTION 

The  heat  of  formation  is  frequently  viewed  as  a  measure  of  the  energy  content  of 
a  compound.  As  such,  it  is  an  important  factor  to  consider  in  designing  new 
energetic  materials  or  evaluating  existing  ones.  For  example,  the  heat  of  formation 
enters  into  the  calculation  of  such  key  explosive  and  propellant  properties  as 
detonation  velocity  [1,2],  detonation  pressure  [1,2]  and  specific  impulse  [3]. 

We  have  recently  developed  a  nonlocal  density  functional  procedure  ior 
computing  gas  phase  heats  of  formation  [4].  The  advantages  of  a  computational 
approach  are  of  course  that  it  can  be  applied  to  suggested  target  molecules  prior  to 
undertaking  syntheses,  and  to  newly-prepared  compounds  even  when  the  available 
amount  is  insufficient  for  laboratory  characterization.  We  have  now  applied  our 
technique  to  a  variety  of  molecules  that  are  of  potential  interest  as  energetic 
materials  [5-8].  These  as  well  as  more  recent  results  will  be  summarized  and 
discussed. 

THEORY 

Our  procedure  [4]  involves  calculating  AE  for  the  formation  of  the  molecule 
from  its  elements  at  0  K  followed  by  conversion  to  AH  at  298  K  by  assuming  ideal 
behavior  and  adding  the  translational,  rotational  and  vibrational  energies. 

Empirical  correction  terms  corresponding  to  the  various  coordination  states  of  the 
carbons,  nitrogens  and  oxygens  are  also  added. 

AE  is  computed  by  a  nonlocal  density  functional  procedure  (Gaussian  94  [9], 

Becke  exchange  and  Perdew  correlation  functionals  [10,11],  6-31G**  basis  set ),  using 
optimized  geometries.  The  vibrational  energy  was  originally  determined  from  the 
normal  mode  frequencies,  the  calculation  of  which  can  be  quite  demanding  in  terms 
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of  computer  resources.  We  have  found,  however,  that  this  energy  can  be  estimated 
accurately  from  stoichiometry-based  relationships  [12];  accordingly  we  now  compute 
the  frequencies  only  when  we  wish  to  confirm  that  the  geometry  corresponds  to  a 
local  minimum  in  the  energy  (indicated  by  the  absence  of  imaginary  frequencies 

7able  I,  our  calculated  heats  of  formation  are  compared  to  the  experimental 
values  for  some 'compounds  relevant  to  the  area  of  energetic  materials.  Additional 
comparisons  have  been  presented  elsewhere  [4].  The  average  absolute  error  m  the 
test  cases  is  3.4  kcal/mole,  indicating  very  satisfactory  agreement. 

The  procedure  that  has  been  described  yields  gas  phase  heats  of  formation,  and 
these  will  be  the  primary  subject  of  this  paper.  For  practical  purposes,  however, 
liquid  and  solid  phase  values  are  often  of  greater  interest;  these  can  be  obtained  if 
the  heats  of  vaporization  and  sublimation  are  known,  using  eqs.  (1)  an  (.-)■ 


1 


? 


a 
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Table  I.  Comparison  of  Some  Calculated  and  Experimental 

Gas  Phase  Heats  of  Formation  _ _ _ 

—  Molecule  Heat  of  formation,  kcal/ mole 


(H3C)2N-N02 


Calculated _ Experimental 

_2.7a  -l.lb 


124.3d 


130e 


125.  ld  124S 


aReference  6.  bReference  45. 

cReference  4.  Reference  8. 

eExtrap  dated  from  the  experimental  values  for  the  dibenzo-  and  monobenzo- 
derivatives,  which  are  142.8  and  136.4  kcal/mole,  respectively,  t 
fy.  t.  Chia  and  H.  E.  Simmons,  J.  Am.  Chem.  Soc.  89,  2638  (1967). 
oExtrapolated  from  the  experimental  values  for  the  dibenzo-  and  monobenzo- 
derivatives,  which  are  132.1  and  128.2  kcal/mole,  respectively.^  _ 
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(1) 

(2) 


AH?98  (liquid)  =  AH^98  (gaseous)  -  AH^'8 
AH?98  (solid)  =  AHp8  (gaseous)  -  AH“98 

We  have  recently  developed  quantitative  representations  of  both  AH^p  an^ 

AH298  terms  of  computed  quantities  related  to  electrostatic  potentials  on 
molecular  surfaces  [14,15];  the  approach  is  one  that  we  have  used  successfully  to 
relate  a  variety  of  liquid,  solid  and  solution  properties  to  quantities  calculated  for 
individual  molecules  [14,16,17].  These  expressions  for  AH^98  an(5  AHj9®  can 
accordingly  be  used  in  conjunction  with  AHP8(gaseous)  to  obtain  liquid  and  solid 
phase  heats  of  formation,  via  eqs.  (1)  and  (2). 

RESULTS  AND  DISCUSSION 

We  will  present  and  discuss  our  computed  gas  phase  heats  of  formation  for 
several  classes  of  compounds  that  are  currently  of  interest  as  potential  energetic 
materials.  Some  of  these  data  have  been  published  elsewhere,  as  will  be  indicated. 
We  will  also  mention  two  solid  phase  heats  of  formation  that  we  have  recently 
obtained,  using  eq.  (2)  and  our  expression  for  AH^g. 

1.  T etraazapentalenes 

The  presence  of  several  linked  nitrogens  in  a  molecule  (nitrogen  catenation)  is 
frequently  associated  with  instability  [18].  It  is  therefore  notable  <hat  certain 
derivatives  of  the  tetraazapentalene  4  show  unexpected  stability  [19-21];  of  particular 
interest  in  the  present  context  is  TACOT,  6,  which  has  a  melting/ decomposition 
temperature  of  378°  C  [22],  We  have  recently  speculated  that  the  surprising  stability 
that  has  been  observed  experimentally  for  the  molecular  framework  4  may  be  due  to 
the  relatively  positive  character  that  we  have  found  for  the  two  triply-coordinated 


In  Table  II  are  listed  our  computed  gas  phase  heats  of  formation  for  the  two 
isomeric  tetraazapentalenes  4  and  5  and  for  three  of  their  derivatives.  None  of  these 
compounds  are  known  at  present,  although  various  other  derivatives  have  been 
prepared  [19-21,23,24],  including  6  [22].  7  [25]  and  8  and  9  [26]  have  been  proposed  as 
target  energetic  materials,  and  synthesis  efforts  have  been  undertaken. 
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Table  II.  Calculated  Gas  Phase  Heats  of  Formation  of  Some  Tetraazapentalenes5 

- - - - - -  '  298 

Molecule 

kcal/mole 


AHf8 

cal /g 


N 


X 


N 


X 


N 


124 


1151 


125 


1158 


253 


1318 


11 

RDX 


no2 

N„ 


45.7  206 

(experimental  valueb) 


0-?N  no2 


aThe  values  for  all  of  the  tetraazapentalenes  except  10  are  from  reference  8. 
^References  27  and  28. _ _ _ _ _ _ _ _ _ 
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The  data  in  Table  II  are  given  in  both  kilocalories /mole  and  calories/ gram;  it  is 
the  latter  that  is  relevant  for  predicting  detonation  performance  [1-3].  In  order  to 
provide  a  basis  for  comparison,  the  experimentally-determined  gas  phase  heat  of 
formation  of  RDX,  11,  a  prominent  military  explosive,  is  also  included  [27,28]. 

Most  of  the  tetraazapentalenes  in  Table  II  have  quite  large  positive  heats  of 
formation,  greater  by  at  least  a  factor  of  five,  in  calories/ gram,  than  that  of 
RDX-  the  hi°h  values  obtained  for  the  unsubstituted  molecules  4  and  5  are 
particularly  striking.  The  relatively  low  AH]  found  for  10  may  indicate  that 


Table  m  Calculated  Gas  Phase  Heats  of  Formation  of  Some  Furazans  and  Tnazoles3- 

Molecule 


798 

AH] 


?98 

AH] 


kcal/ mole 


cal/ 


40 


202 


330 


903 


14 

DNAF 


N  N  O  N  N 
O  O 


169 


621 


15 

NTO 


-5.7  -44 


b  References  27  and  28. 


k/i 
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framework  nitrogens  are  more  effective  than  nitro  groups  in  producing  a  large  heat 
of  formation.  The  diazide  8  and  the  ditetrazole  9  are  tautomers,  which  are  expected 
to  be  related  through  a  tautomeric  equilibrium,  8  f=^  9.  Our  results  show  the 
diazide  to  be  the  more  stable,  by  about  16  kcal/ mole.  This  could  be  a  matter  for 
concern,  from  the  standpoint  of  an  energetic  material,  because  of  the  tendency  of 
many  azides  toward  facile  decomposition  [29]. 

2.  Furazans  and  Triazoles 

Like  the  tetraazapentalenes,  the  furazans  and  triazoles  in  Table  III  fit  into  the 
category  of  high-nitrogen  compounds.  However  these  have  actually  been 
synthesized:  12  [30,31],  NOTO,  13  [32],  DNAF,  14  [33],  NTO,  15  [34,35],  ADNT,  16  [36]. 
NOTO  has  the  interesting  feature  of  being  a  liquid  at  room  temperature  [32].  Table 
III  shows  that  all  of  these  compounds  except  NTO  have  positive  heats  of  formation 
which  are  considerably  larger,  on  a  calories/ gram  basis,  than  that  of  RDX  but  much 
smaller  than  those  of  most  of  the  tetraazapentalenes  in  Table  H. 


3.  Nitra mines  and  Difluoramines 


The  difluoramino  group,  — NF2,  continues  to  be  of  interest  as  an  ingredient  of 
energetic  materials,  particularly  propellants  [2,3,37,38].  A  judicious  combination  of 
nitro  and  difluoramino  substituents  offers  the  possibility  or  enhanced  propellant 
performance  due  to  an  increased  number  of  moles  of  gaseous  combustion  products 
per  gram  of  material  (provided  that  some  hydrogens  are  present)  [3,6,39]. 

In  conjunction  with  recent  computational  studies  of  various  existing  or 
proposed  nitramine  and  difluoramine  systems  [6,40],  we  have  calculated  the  gas 
phase  heats  of  formation  of  the  compounds  in  Table  IV.  While  most  of  them  are 
positive,  only  18  has  a  gas  phase  value  (in  calories/ gram)  greater  than  that  of  RDX. 

A  consistent  pattern  in  Table  IV  is  that  AH*[  is  lower  (less  positive)  for  the 
difluoramine  than  for  its  nitramine  analogue.  (Leroy  et  al  have  observed  the  same 
to  be  true  when  the  -NF2  and  -NCb  are  attached  to  carbons  [41],  as  we  do  also  for  the 
pair  20,  22.)  However  the  difluoramine  may  still  have  a  higher  specific  impulse  (a 
measure  of  propellant  thrust  [3,22,39])  due  to  producing  more  moles  of  gaseous 
products  per  gram  upon  combustion.  Thus  the  specific  impulse  of  2j  is  predicted  to 
be  higher  than  that  of  RDX,  11,  despite  having  a  lower  estimated  heat  of  formation 
[3]. 


23 


NF, 
I  " 


02N  no2 


O' 
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Table  IV.  Calculated  Gas  Phase  Heats  of  Formation  of  Some 

Nitramines  and  Difluoramines3 

Molecule 

AHf8 

kcal/mole 

298 

AHf 
cal/ g 

1 

(H3C)2N-N02 

-2.7 

-30 

17 

(H3C)2N-NF2 

-20 

-206 

18 

02N — N^N—  N02 

44 

297 

19 

02N — NF2 

29 

186 

20 

TNAZ 

Oi\/\ 

V  N— N02 

o2n  ^ 

31 

160 

21 

02N 

X  N— NFZ 

01N/A^ 

16 

79 

22 

OiK/\ 

x  N—  N02 

7.4 

38 

11 

RDX 

45.7  206 
(experimental  valueb) 

a  All  of  these  calculated  results  except  that  for  22 

are  from  reference  6. 

b  References  27  and  28. 

- - - 

4.  Some  C.  N.  0-  F  Compounds. 

We  have  investigated  computationally  a  group  of  19  unsaturated  molecules  that 
contain  only  the  elements  C,  N,  O  and  F  (Table  V)  [7],  Their  high  N/C  and  O/C 
ratios  suggest  that  these  are  potential  energetic  molecules  or  precursors.  To  the  best 
of  our  knowledge,  most  of  these  molecules  are  presently  unknown,  an  exception 

being  26  [42,43]. 

The  calculated  gas  phase  heats  of  formation  are  given  in  Table  V.  These  were 
obtained  by  a  slightly  modified  version  of  the  method  described  earlier,  in  that  ab 
initio  HF /6-31G*  optimized  geometries  and  vibrational  frequencies  were  used  and 
AE  was  computed  with  the  density  functional  program  DeMon  [44].  Structures  26A, 
dinitrosoacetylene,  and  26B,  which  is  more  accurately  designated  as  dinitrile-di-N- 
oxide,  were  both  found  to  correspond  to  energy  minima  at  the  ab  initio  HF/6-31G* 
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Table  V.  Calculated  Gas  Phase  Heats  of  Formation  of  Some  C,  N,  Q, 
-  -  ~  AH?98 


Molecule _ 

o2n-o=c-no2 

o,n-<=c-no 


CN 

(0,N),C — C — O 

o2n— n=c=n— no, 

(02N)2C  =t_=C=0 
(0,N),C=C=C=N-N0, 
0,N— N=C=C=N— NO? 


0,N-N= 
0,N— N- 


0,N  NO, 

V 

/_\_ 

O,^  NO, 


:=n-no, 


kcal/mole 

89 

100 


F  Compounds3 
AHf8 
cal/ g 


26A 

ON— C=C— NO 

130 

1550 

26B 

O 

£ 

III 

Y 

V 

Z 

O 

73 

868 

27 

0,N — C^C  NF, 

74 

606 

28 

ON — C=C- — NF, 

95 

901 

29 

N3-G=C— N02 

144 

1287 

30 

OCN — C^C  NO, 

54 

486 

NO, 

80 

426 

31 

(0,N),C=C^ 

CN 

32 

0,N 

^C=C(NO,),_ 

108 

572 

N-NO, 


0,N  NO, 

C(N02)2 

c 

/\ 

/C==CS 

0,N  NO, 


(continued) 
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Table  V.  Calculated  Gas  Phase  Heats  of  Formation  of  Some  C,  N,  O,  F  Compoundsa 

( continued) _ 


Molecule 


AH?98 

kcal/mole 


AHf8 
cal/ g 


O 


42 


(02N)2Cx  7C(N02)2 

c 

II 

o 


35 


11 


RDX 


45.7  206 

(experimental  valueb) 


aThe  results  for  24  -  30,  33-35,  37  and  38  are  in  reference  7. 
b  References  27  and  28.  _ _ _ 


level  [7].  Table  V  shows  that  26B  is  the  more  stable,  which  is  consistent  with 
experimental  evidence  [42,43]. 

Most  of  these  compounds  are  predicted  to  have  rather  large  positive  heats  of 
formation  (calories /gram).  One  of  the  highest,  not  surprisingly,  is  the  azide  29.  The 
pairs  24,  27  and  25,  28  show  again  that  replacing  -NO2  by  -NF2  lowers  AHf98,  as 
observed  above  when  these  groups  were  on  nitrogens.  ^ 

Particularly  striking,  especially  on  a  calories/ gram  basis,  is  the  increase  in  AHr 
that  occurs  when  the  nitro  group  is  replaced  by  the  nitroso.  This  can  be  seen  in  24, 
25  and  26 A  and  again  in  27  and  28.  (26B  is  more  nitrile-N-oxide  than  nitroso  in 
character,  as  mentioned  above.)  This  effect  has  also  been  observed  experimentally; 
the  heat  of  formation  of  43  is  46.44  kcal/mole  (322.5  cal/g)  vs.  13.9  kcal/mole  (79.0 
cal/g)  for  44  [45].  In  view  of  these  findings,  it  is  interesting  to  note  evidence 


NO 

I 

,N. 

1 

NO 

43 


N02 

'N 

1 

no2 

44 


0,N 


NO 

rS 


"NO-, 


45 


N02 
I  z 

A 

N'  N 

0,N  N02 


11,  RDX 


suggesting  the  possibility  of  another  desirable  consequence  of  substituting  -NO  for 
-NO2;  an  experimental  study  showed  that  the  impact  sensitivity  of  45  is  significantly 
less  than  that  of  RDX,  11,  while  the  explosive  power  (measured  by  dent  depth)  is 
essentially  the  same  [46]. 

A  further  comparison  of  interest  is  between  =C(N02)2  compounds  and  their 
=N(NO?)  analogues.  The  latter  have  the  greater  heats  of  formation,  as  can  be  seen 
from  the  pairs  36,  37  and  40, 41. 
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5.  Solid  Phase  Heats  of  Formation 

As  was  pointed  out  earlier  in  this  paper,  we  now  have  the  means  for  reliably 
estimating  heats  of  vaporization  and  sublimation  [14,15].  These  permit  us  to 
convert  our  computed  gas  phase  heats  of  formation  to  liquid  and  solid  phase  values, 
using  eqs.  (1)  and  (2),  respectively.  We  have  thus  far  applied  this  procedure  to  two 
compounds,  1,3,3-trinitroazetidine  (20,  TNAZ)  and  DNAF  (14).  ^ 

For  TNAZ,  we  find  AHgg  =  22  kcal/mole;  when  combined  with  AH'  (gaseous) 

=  31  kcal/mole  (Table  IV)  via  eq.  (2),  the  result  is  AH]9  8(solid)  =  9  kcal/mole.  This  is 
in  excellent  agreement  with  the  experimental  value,  8.7  kcal/mole  [38].  For  DNAF, 
we  obtain  AH?ub  =  33  kcal/mole.  Using  AH']-  (gaseous)  =  169  kcal/mole  (Table  HI), 

9QQ 

eq.  (2)  gives  AH]  (solid)  =  137  kcal/ mole. 

CONCLUSIONS 

We  have  presented  and  discussed  the  results  of  density  functional  calculations  of 
gas  phase  heats  of  formation  for  several  categories  of  energetic  compounds.  We 
have  also  demonstrated  the  feasibility  of  converting  these  to  solid  phase  values. 

For  the  systems  investigated,  some  key  observations  are: 

(1)  The  tetraazapentalenes  have  particularly  large  positive  heats  of  formation 

(calories/  gram). 

(2)  Replacement  of  a  nitro  by  a  difluoramino  group,  whether  on  a  carbon  or  on  a 
nitrogen,  lowers  the  heat  of  formation. 

(3)  Replacement  of  a  nitro  by  a  nitroso  group  increases  the  heat  of  formation. 

(4)  Replacement  of  =C(N02)2  by  =N— NO2  increases  the  heat  of  formation. 
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